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ABSTRACT 


This  report  describes  the  work  accomplished  on  the  retrieval  of  strato¬ 
spheric  temperature  and  density  profiles  from  remote  earth  limb-view  infrared 
radiation  measurements .  Temperature  or  density  can  be  retrieved  in  a  routine 
manner.  A  model  based  on  two  rotational  lines  of  CO2  was  developed  to  retrieve 
simultaneously  both  temperature  and  density.  This  model  is  running  on  a  CRAT-1 
Computer  but  needs  additional  verification.  An  investigation  of  alternative 
methods  to  calculate  infrared  transmittance  has  showed  that  an  interpolation 
scheme  does  work  for  limb-view  geometries  with  a  subsequent  potential  savings  in 
computer  time  of  one  thousand  or  better.  Though  the  concept  has  been  proved,  the 
ultimate  accuracy  of  this  method  has  still  to  be  evaluated. 


I.  INTRODUCTION 


This  report  describes  the  work  that  was  done  on  the  project  "Atmospheric 
Temperature  profiles  from  Earth  Limb  Infrared  Measurements"  for  the  Air  force 
Office  of  Scientific  Research,  grant  no.  AFOSR  -  82  -  0153.  Numerical  methods 
were  evolved  to  utilize  remote  sensing  of  the  earth  limb  by  infrared  radiation 
measurements  in  order  to  determine  atmospheric  temperature  and  constituent 
profiles.  These  methods  are  primarily  applicable  to  stratospheric  altitudes  and 
for  IR  measurements  such  as  the  AFGL  Cryogenic  Infrared  Radiance  Instrumentation 
for  Shuttle  (CIRRIS).  The  work  covered  two  principal  areas:  1)  An  inversion 
scheme  utilizing  two  rotational  lines  of  CO2  to  retrieve  simultaneously  tempera¬ 
ture  and  density  profiles  between  about  10  and  65  kms;  2)  A  method  to  calculate 
atmospheric  transmittance.  This  latter  task  was  subdivided  into  the  jobs  of  dev¬ 
eloping  a  very  accurate  scheme  to  calculate  transmittance  of  a  single  rotational 
line  of  C02  and  then  proving  the  feasability  of  an  alternate  interpolation 
scheme  to  calculate  transmittance  that  would  be  extremely  fast  and  take  less 
computer  time.  The  details  of  the  inversion  algorithm  are  discussed  in  the  next 
section.  Conclusions  and  suggestions  for  further  work  are  given  in  the  third 
section.  Theory  covering  radiation  transport  is  presented  in  the  appendix. 

II.  TECHNICAL  DISCUSSION 

Prior  to  the  start  of  this  work  the  author  had  developed  two  computer  codes 
which  would  take  a  set  of  earth  limb  IR  radiance  measurements  and  calculate  a 
vertical  temperature  or  density  profile  for  the  radiating  species.  (Figure  1 
illustrates  the  limb  viewing  geometry.)  The  method  was  essentially  one  of 
choosing  some  initial  profile  (temperature  or  density)  and  calculating  the  cor¬ 
responding  limb  radiance.  Using  the  difference  between  the  measured  and  calcu¬ 
lated  radiances  as  a  modification  to  the  initial  profile,  a  second  and  hopefully 
better  profile  could  be  obtained!  By  successive  iterations  the  first  guess 

+  Chatline  M.  T.,  A  General  Relaxation  Method  for  Inverse  Solution  of  the  Full 
Radiative  Transfer  Equation,  J.  Atmos.  Scl.,  29*  pp.  741,  1972 
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Geometry  for  Lnfrared  Limb  Radiance  Inversion  Calculations 


profile  could  be  reduced  to  a  profile  which  radiated  with  the  same  altitude 
intensities  as  those  measured.  Figure  2  illustrates  the  convergence  process 
suggested.  Thus  far  we  have  not  used  real  radiation  data  since  none  existed  of 
sufficiently  fine  resolution.  Our  "measured"  data  was  actually  calculated  from 
a  reference  atmosphere.  For  temperature  retrieval  the  density  had  to  be  known. 
For  density  retrieval  the  temperature  had  to  be  given.  The  theory  of  radiation 
transfer  as  applied  to  this  remote  sensing  is  given  in  Appendix  A. 

Two  Line  Model 

This  work  of  temperature  and  density  retrieval  could  be  directly  applicable 
to  the  AFGL  Cryogenic  Infrared  Radiance  Instrumentation  for  Shuttle  (CIRRIS) 
data  because  the  wavelength  resolutions  are  about  the  same.  However,  with  CIRRIS 
data  it  will  not  be  possible  to  specify  either  a  temperature  or  density  for 
retrieval  of  the  other,  tfhat  is  needed  is  a  method  to  retrieve  both  temperature 
and  density  simultaneously  in  a  3elf  consi3tant  manner.  The  problem  was  attacked 
by  first  utilizing  the  existing  computer  codes  in  a  push-pull  arrangement.  A 
first  guess  for  both  temperature  and  density  profiles  was  used  to  calculate  the 
radiances  for  one  rotational  line.  The  rotational  lines  of  the  CC^  15um  band 
were  used  for  this  work  since  the  transmittance  properties  vary  from  optically 
thin  to  very  optically  opaque.  (The  whole  line  was  used  initially  but  the  trans¬ 
mittance  subroutine  has  now  been  modified  to  handle  any  instrument  slit  width.) 
The  difference  between  the  calculated  and  measured  radiations  was  used  to  modify 
the  initial  temperature  profile.  This  new  temperature  profile  and  the  initial 
density  guess  were  then  used  to  calculate  radiation  from  the  second  line.  The 
difference  between  calculated  and  measured  radiations  from  this  second  line  was 
used  to  modify  the  density.  Both  the  new  density  and  temperature  profiles  were 
now  used  to  calculate  a  new  value  of  radiation  intensity  for  the  first  line. 
Successive  radiation  calculations  back  and  forth  between  the  lines,  together 
with  the  successive  modifications  of  the  temperature  and  density  profiles  was 
supposed  to  converge  simultaneously  the  temperature  and  density  profiles  to 


Convergence  of  Chatline* s  Relaxation  Method  for  CO-  Density  Retrieval 


unique  values  compa table  with  the  radiation  measurements.  Unfortunately  the 
method  did  not  perform  as  anticipated.  Temperature  and  density  profiles  would 
indeed  converge  from  their  original  values  to  within  about  10$  of  the  true 
values,  but  then  for  some  reason  the  temperature  and  density  profiles  diverged 
away  from  the  expected  values.  This  behavior  was  very  unsettling  and  completely 
unexpected  after  the  absolute  convergence  of  our  earlier  results.  Because  of  the 
large  computer  costs  to  run  this  double  line  code  we  were  not  able  to  explore  or 
to  understand  the  details  of  this  convergence  behavior,  life  thought  that  the 
trouble  was  caused  by  inaccuracies  in  the  transmittance  calculation.  Accordingly 
a  major  effort  was  undertaken  to  improve  the  accuracy  of  the  integrated  trans¬ 
mittance.  Those  results  are  discussed  in  the  next  section.  At  about  this  time  a 
new  two  line  retrieval  scheme  was  also  devised.  This  scheme  is  based  on  the 
Taylor’s  series  evaluation  of  a  function.  If  a  quantity  F  is  a  function  of  two 
independent  variables  x  and  y,  then  the  Taylor' 3  series  expansion  is 
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where  h  -  x-xQ  and  k  -  y-y^.  if  h  and  k  are  sufficiently  small  then  the  first 
two  terms  of  the  series  form  a  reasonable  approximation.  Since  the.  limb  radi¬ 


ation  values  are  functions  of  both  the  temperature  and  density,  the  above 

/ 

expression  can  be  used  to  write 


F(T,  D) 
G(T,  D) 


where  F  and  G  are  the  measured  radiances  of  twb  C0o  lines,  T  and  D,  the  first 

c.  0  0 

estimate  of  the  temperature  and  density  values  which  give  the  trial  intensities 
F(Te,  Dq)  and  G(T0,  Dq),  and  T  and  0  are  the  new  temperature  and  density  values 
which  are  a  linear  extrapol  Mon  for  le  measured  intensities  F(T,  D)  and 
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G(T,  D).  The  matrix  of  partial  derivatives  is  called  the  Jacobian  matrix.  Rear¬ 
ranging  and  multiplying  by  [j]”1  we  get  the  iteration  formula. 
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where  INTI  and  INT2  are  the  measured  IR  radiation  values.  The  advantage  of  this 
Jacobian  matrix  method  is  that  the  convergence  proceeds  along  the  path  of 
maximum  descent  and  does  30  simultaneously  for  both  temperature  and  density.  The 
method  has  been  implimented  as  shown  in  the  flow  chart  of  Figure  3.  This  method 
was  coded  up  for  the  NCAR  CRAY-1  Computer,  but  we  ran  out  of  computer  funds 
before  the  debugging  was  completed.  Discussions  with  the  mathematicians  assure 
us  that  tie  scheme  will  converge  provided  our  initial  guesses  for  TQ  and  Dq  are 
not  too  far  away  from  the  desired  values. 

Transmittance  Calculations 

A  calculation  of  the  atmospheric  transmittance  is  necessary  whenever  a 
value  of  the  limb  radiation  is  needed.  Infrared  radiation  in  the  atmosphere  i3 
caused  by  rotational  transitions  in  molecules  which  have  a  magnetic  moment. 
Because  the  molecule  has  translational  kinetic  energy  and  is  buffeted  by  other 
ir  ‘ ccules  while  in  a  rotational  transition,  the  absorption  -  emission  line  shape 
i.e.  the  variation  of  the  radiation  intensity  versus  frequency,  changes  de¬ 
pending  upon  its  position  in  the  atmosphere.  3elow  35  km  the  line  is  collision 
broadened,  from  40  to  150  km  the  line  is  doppler  broadened.  Since  these  widths 
are  different  one  must  account  for  the  variable  line  shape  when  calculating 
transmittance.  The  only  method  that  currently  exists  to  handle  the  variation  in 
line  width  is  to  calculate  transmittances  at  fixed  frequency  values  distributed 
across  the  line  width.  By  watching  how  these  individual  transmittances  vary 
through  the  absorbing  medium  we  can  then  determine  the  average  transmittance 
characteristic  of  the  entire  line.  Not  only  does  the  line  width  vary  with 


altitude  due  to  collisions  and  motion  but  the  amplitude  of  each  line  relative  to 


Ill 


the  other  lines  changes  with  the  change  of  atmospheric  temperature.  This  is 
reflected  in  the  changes  in  the  shapes  of  the  P  and  R  branches  as  a  function  of 
temperature . 

A  Voigt  profile,  which  is  a  convolution  of  the  Doppler  and  Lorentz 
(collision)  profiles,  was  used  to  describe  the  smoothly  changing  line  shape.  An 
emperical  formula  based  on  the  ratio  of  the  Lorentz  width  to  doppler  width  was 
devised  to  specify  the  effective  line  width  for  all  altitudes.  Forty  five  points 
from  the  line  center  to  half  way  between  lines  was  considered  adequate  frequency 
resolution.  Fifteen  points  were  used  to  specify  the  core  of  the  line,  10  more 
points  through  the  wings,  and  up  to  20  more  to  reach  the  frequency  midway 
between  lines.  Forty  five  points  are  probably  more  points  than  are  absolutely 
necessary  but  we  wanted  to  be  on  the  conservative  side.  As  the  rotational  line 
broadens  and  flattens  out  less  points  are  required  to  accurately  specify  the 
shape.  Figure  4  is  a  plot  of  the  transmittance  along  a  limb  view  line  of  sight 
from  outside  the  atmosphere  into  the  tangent  altitude  for  several  such  optical 
paths.  These  curves  were  all  calculated  for  the  J  ■  26  rotational  line  of  the 
15um  C02  P  branch.  Note  the  complete  absorption  at  the  center  of  the  line  and 
the  broadening  with  decreasing  altiitude.  Also  shown  are  the  spacings  of  the 
frequency  points  through  the  wing  of  the  line.  The  smooth  shape  and  the  coverage 
afforded  by  the  frequency  subdivisions  should  be  accurate  enough  to  handle  all 
the  forseeable  IR  data. 

Any  infrared  sensor  has  a  vertical  resolution  which  depends  upon  the  size 
of  the  telescope.  For  the  SPIRE  instrument  this  vertical  resolution  was  an  8  km 
footprint  at  the  tangent  heights.  For  the  CIRRIS  instrument  and  for  probably  all 
future  data  a  3  km  footprint  size  /as  deemed  more  satisfactory.  It  was  found 
that  radiation  calculations  depended  critically  on  the  way  the  density  was 
averaged  over  the  3  km  interval.  Vfe  finally  determined  that  transmittances 
should  be  calculated  for  each  1  km  thick  optical  path  and  that  a  plain  average 
for  the  three  1  km  thick  transmittances  would  be  used  for  the  radiation  calcula- 


Transmittance  at  the  Indicated  Tangent  Altitudes  Showing  the  Frequency  Distribution 
and  Line  Shape 


tion.  (See  Figure  1.)  Within  each  3  km  footprint  the  density  was  allowed  to  have 
an  exponential  variation  with  height.  Calculations  revealed  that  radiations  from 
3  km  thick  slabs  using  the  average  transmittance  were  almost  identical  to  the 
sum  of  the  radiations  from  3  of  the  1  km  thick  slabs.  This  must  mean  that  the 
integrated  transmittance  over  the  entire  line  has  a  much  less  extreme  variation 
with  density  than  the  individual  monochromatic  transmittances.  An  average  temp¬ 
erature  of  the  three  1 -km  slabs  was  found  to  be  adequate  for  the  calculations. 
The  iteration  scheme  returns  a  single  temperature  and  density  value  for  each 
3  km  slab  of  atmosphere.  Within  each  3  lcm  interval  along  the  line  of  sight  the 
change  in  transmittance,  ^Tr,  was  found  for  each  individual  frequency.  The  total 

monochromatic  transmittance  at  frequency  v .  from  level  i  to  the  top  of  the 

J 

atmosphere  could  then  be  written. 

Tri+i  (v.)  -  ATr.U.)  Tr.(v.) 

Since  the  line  width  changes  with  altitude,  the  previous  values  of  transmittance 
Tr^Cv^)  had  to  be  interpolated  to  the  new  frequencies  v^.  The  integrated  trans¬ 
mittance  was  then  found  by  integrating  Tri+1  (v^)  over  frequency.  The  integrated 
transmittance  was  given  numerically  by 

45 

Tri.t  ‘  5-i  Tri*i  (vj)  avi 

The  variation  of  transmittance  for  the  J  *  26  rotational  line  of  the  C02  15pm  P 
branch  is  plotted  for  eight  different  tangent  altitudes  in  Figure  5«  Obviously 
the  total  transmittance  increases  for  increased  depth  into  the  atmosphere.  The 
radiation  leaving  the  atmosphere  along  the  i^*1  tangent  line  of  sight  is  found  by 
the  formula 

ft 

-  /  B(T)  dTr  (T,  C02,  s) 

where  B  is  the  Planck  Function  and  the  transmittance  Tr  is  a  function  of  temper¬ 
ature,  density  and  obviously  position  along  the  line  of  sight.  According  to 
this  expression  the  maximum  radiation  comes  from  the  altitudes  where  the  trans- 


mittance  has  the  greatest  change.  Figure  6  is  a  plot  of  transmittance  for  a 
number  of  different  rotational  lines  and  different  tangent  altitudes.  It  can  be 
seen  that  the  radiation  at  different  lines  comes  from  different  altitudes.  For 
the  simultaneous  temperature  and  density  retrieval,  we  chose  the  J  *  26  line 
which  has  good  variation  between  1 1  and  65  kms  for  one  of  the  lines.  The  other 
line  was  a  composite  of  the  J  *  2  and  J  *  31  lines.  Between  11  -  17  km  the 
J  *  31  line  gave  adequate  change  in  transmittance,  while  the  J  ■  2  line  was 
satisfactory  between  20  and  65  km. 

The  method  for  calculating  transmittances  a3  briefly  described  above  is 
accurate  but  takes  an  excessive  amount  of  computer  time.  For  real  data  sets 
such  as  will  be  obtained  from  the  CIRRIS  experiment  a  much  faster  and  cheaper 
scheme  must  be  found  to  calculate  transmittance.  We  believe  such  a  scheme  is 
afforded  by  the  McMillin*  technique. 

McMillin  Calculation  of  Transmittance 

McMillin’ s  method  was  designed  and  has  been  successfully  used  for  inversion 
of  data  from  vertically  downward  looking  instruments  which  have  a  band  width  of 
many  rotational  lines  and  the  band  includes  radiation  from  several  different 
species.  For  the  limb  viewing  geometry  we  desire  a  calculation  of  transmittance 
over  a  bandwidth  of  only  one  rotation  lice  and  considering  only  one  species. 
The  transmittance  is  calculated  in  steps  along  the  line  of  sight  path  using  the 
expression 

Tr.  -  ^Tr.  Tr.  , 
i  i  i-l 

where  the  transmittance  is  given  by  the  expansion 

12 

i,ri  ■ 

*  McMillin,  L.  M.,  H.  F.  Fleming  and  M.  L.  Hill,  Atmospheric  transmittance  of  an 
absorbing  gas.  3:  A  computationally  fast  and  accurate  tranmittance  model 
for  absorbing  gases  with  variable  mixing  ratios,  Appl.  Optics,  18, 
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Figure  6.  Variation  in  Transmittances  for  the  CC^  Rotational  Lines  and  Tangent 
Altitudes  Shown 
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The  hat  ~  over  T  and  p  refers  to  reference  atmosphere  values.  In  the  calculation 

/v 

a  reference  atmosphere  must  be  chosen  with  temperature  profile  T  and  pressure 

A 

profile  p.  The  total  column  density  for  the  reference  atmosphere  is  larger  than 
for  any  other  atmosphere  for  which  the  transmittance  is  to  be  calculated.  The 
independent  variable  in  these  expressions  is  column  density  along  the  path 

T.  oofvn(Z)da 

which  is  also  called  the  optical  depth  in  (l)  of  the  appendix.  The  various 
starred  quantities  along  the  path  are  weighted  averages  defined  by 


r 


Ap  -  p  -  p 

In  order  to  use  the  transmittance  as  given  by  the  series  expansion  it  is 
necessary  to  have  ATr^  small.  HcMillin  suggested  that  100  steps  of  about  equal 
transmittance  difference  should  give  good  accuracy.  To  generate  the  coefficient 
j ,  a  number  of  model  atmospheres  must  be  chosen  and  used  as  the  basis  atmo¬ 
spheres.  The  basis  must  contain  more  than  12  different  models  and  the  models 
should  be  chosen  to  cover  as  wide  a  range  of  parameters  as  is  expected  in  the 
transmittance  calculations.  These  basis  models  are  used  to  generate  transmit¬ 
tance  values  at  each  of  the  reference  column  densities  t^.  The  C^'s  are  found 
which  give  the  best  mean  square  fit  to  these  transmittance  of  the  basis  atmo¬ 
spheres.  The  exception  to  this  mean  square  fit  is  C , .  which  is  chosen  so  that 


Tri-1 

Our  test  of  the  McMillin  calculation  used  the  J  *  10  rotational  line  of  the 
P  branch  of  the  15um  band  of  COg.  The  model  atmospheres  for  the  basis  atmosphere 
were  generated  mathematically  using  the  hydrostatic  equation,  constant  00^ 
mixing  ratio  (3.16x10“^),  and  various  assumed  temperature  profiles.  Some  of 
these  profiles  are  shown  in  Figure  7.  All  the  temperature  profiles  fell  within 
♦10°K  of  the  solid  line  profile  of  Figure  7.  A  total  of  nine  different  profiles 
were  used  to  generate  nine  different  atmospheres.  Sach  of  these  model  atmo¬ 
spheres  was  used  to  generate  transmittance  curves,  such  as  shown  in  Figure  4, 
for  each  tangent  height.  Actually  there  were  17  tangent  heights  of  3  km  slabs 
between  12  km  and  60  km.  The  reference  model  atmosphere  was  chosen  from  the 
basis  atmospheres  to  be  that  atmosphere  with  the  maximum  column  density.  The 
model  atmosphere  to  the  extreme  right  of  Figure  7  was  the  reference  atmosphere 
for  these  calculations.  Using  this  atmosphere  and  the  smallest  tangent  height, 
a  calculation  of  transmittance  vs.  column  density  was  made.  The  reference  column 
density  was  found  by  dividing  the  total  transmittance  into  100  equal  multlpli- 
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Figure  7.  Temperature  Profiles  of  Basis  Atmospheres  used  for  McMillin  Transmittance 
Calculations. 


cative  steps  i.e.  malting  the  ATr^'s  in  the  formula 

-  Tri-t  aT^ 

all  equal  with  equal  to  the  maximum  transmittance  calculated.  The 

positions  corresponding  to  this  set  of  reference  Tr^'s  also  specify  a  corre¬ 
sponding  set  of  t^'s.  (The  imply  a  set  of  reference  pressure  and  temp- 

A  A  A  A 

erature  T^.)  The  Tr^'s  can  he  used  immediately  to  generate  *  Tr^/Tr^  1 .  To 

find  the  other  C.  . 's  we  calculate  the  transmittance,  the  difference  between 

J 

A 

temperatures  a^  *  T^-T^,  and  the  pressure  difference  Ap^  *  Pi~^i*  at  Ti'3 
for  each  of  the  other  atmospheres  and  at  all  tangent  heights.  This  gives  us  a 
set  of  A^r^  at  each  which  then  allows  a  mean  square  fit  to  evaluate  the 


s. 


One  problem  with  this  fitting  technique  occurs  because  the  different  tan¬ 
gent  height  calculations  for  the  same  model  atmosphere  differ  only  significantly 
near  the  tangent  point.  This  causes  an  abnormal  weighting  of  the  least  square 
fit  of  ATr^.  To  alleviate  this  problem,  if  the  set  of  ATr^'s  at  some  contain 
a  number  of  ATr^'s  which  are  identical,  only  two  of  the  identical  values  are 
used  in  the  fit. 


With  the  C^'s  calculated,  it  was  now  time  to  evaluate  the  McMillin 
expansion.  Accordingly  a  test  model  atmosphere,  shown  as  the  dotted  profile  in 
Figure  8,  was  devised  to  lie  within  the  limits  of  the  basis  atmospheres  but  not 
equal  to  any  of  them.  The  exact  tranmittance  was  calculated  with  the  43  point 
routine  described  earler.  The  HcMillin  transmittances  were  found  by 
substituting  the  values  of  temperature  from  Figure  8  at  the  reference  t^'s  into 
aT^  and  the  pressure  into  ap^«  The  corresponding  values  of  x^  were  then  found 
at  each  of  the  100  positions.  Finally  the  expansion 


ATr, 


was  evaluated  and  the  transmittances  calculated  by 


Figure  8,  Model  Atmosphere  to  Test  McMillin  Expansion  for  Transmittance. 


Table  I  shows  a  comparison  of  the  exact  transmittance  and  the  transmittance 
calculated  from  the  McMillin  formulas.  Also  shown  is  the  absolute  error.  Mote 
the  absolute  errors  are  very  small,  but  they  are  also  rather  difficult  to 
interpret.  One  possible  way  of  determining  their  significance  is  to  compare 
them  with  the  change  in  transmittance  caused  by  a  change  in  temperature  (the 
sensitivity  of  transmittance  to  temperature).  This  estimate  of  sensitivity  is 
also  shown  at  steps  25  and  40  along  the  line  of  3ight. 

The  calculations  of  transmittance  for  the  test  model  atmosphere  was  found 
to  give  poor  results  beyond  step  40.  This  was  caused  by  too  few  different  atmo¬ 
spheres,  in  our  case  only  9,  used  in  the  determination  of  the  C^'s  at  large 
t^'s.  At  least  12  different  ATr^'a  at  every  are  needed  to  make  a  good  mean 
square  fit.  idith  9  model  atmospheres  and  17  tangent  heights  this  is  easy  for  the 
case  of  small  t^'s  but  not  large  t^'s  since  only  one  tangent  height  produces  the 
largest  A  basis  with  perhaps  20  model  atmospheres  should  correct  this 

problem. 

Note  that  the  McMillin  formulation  that  has  been  discussed  above  attempts 
to  fit  the  transmittance  calculation  at  all  tangent  heights  with  one  set  of 
coefficients  C^.  This  was  only  partially  successful.  At  the  higher  tangent 
heights  the  column  densities  are  such  that  even  the  first  step  of  transmittance 
is  not  seen  along  the  entire  path.  Perhaps  the  calculations  for  different 
tangent  heights  will  have  to  be  handled  with  different  sets  of  coefficients 
based  on  different  column  density  steps  t^.  It  does  look  as  if  a  group  of 
tangent  heights  can  be  handled  with  the  same  set  of  column  density  steps  and 
coefficients  Sven  if  a  set  of  coefficients  were  needed  for  each  tangent 

altitude,  the  computer  storage  requirements  would  be  small  enough  to  fit  on  a 
VAX  type  computer. 


TABLE  I 


Comparison  of  Exact  and  McMillin  Calculation 
Of  Transmittance  (Tri)  Tangent  Height  —  13  km. 


Step 

(i) 

(Tr  ) 
Exact 

(Tr.) 

McMillin 

Aba 

Error 

Aba 

Sensitivity 
(per  °K) 

5 

0.96849 

0.96843 

0.00006 

10 

0.93732 

0.93776 

0.00044 

15 

0.90870 

0.90803 

0.00067 

20 

0.87969 

0.87901 

0.00068 

.00017 

25 

0.85506 

0.85397 

0.00109 

30 

0.82840 

0.82803 

0.00037 

35 

0.79827 

0.79809 

0.00018 

40 


0.74487 


0.74781 


0.00594 
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III.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  NORK 


A  model  of  infrared  radiation  from  CC^  has  been  developed  that  is  appli¬ 
cable  to  earth  limb  radiation  from  the  stratosphere.  The  kernel  of  this  model  is 
a  numerical  code  to  calculate  the  integrated  transmittance  of  a  single  rota¬ 
tional  line.  (However,  a  modification  has  been  made  so  that  any  fraction  of  the 
rotational  line  can  also  be  handled.)  Using  a  non-linear  iteration  process  it 
has  been  shown  that  temperature  or  CO2  density  profiles  in  the  stratosphere  can 
be  found  from  IR  limb  radiation  measurements.  A  numerical  code  has  been 
developed  to  utilize  two  rotational  lines  of  CC^  to  retrieve  both  temperature 
and  CC>2  density  simultaneously.  This  code  is  running  but  need3  verification.  If 
CO2  is  well  mixed  in  the  atmosphere  a  simultanious  retrieval  of  temperature  and 
density  will  yield  the  atmospheric  pressure  profile  and  allow  numerous  densities 
of  other  IR  active  species  to  be  deduced. 

A  computationally  fast  method  for  evaluation  of  transmittance  has  also  been 
studied.  This  interpolation  scheme,  first  proposed  by  McMillin  for  vertical 
looking  sensors,  has  been  reformulated  for  limb  viewing  and  run  with  a  submini¬ 
mum  number  of  basis  atmosphere  models  to  prove  out  the  concept.  Results  were 
very  encouraging.  An  evaluation  of  this  technique  with  12  —  20  model  atmospheres 
should  be  made  to  see  if  one  set  of  expansion  coefficients,  which  only  Lave  to 
be  calculated  once  for  each  sensor  and  each  radiation  band,  can  be  used  for  the 
entire  atmosphere  or  the  other  extreme  of  a  set  of  coefficients  for  each  tangent 
altitude.  The  ultimate  accuracy  that  can  be  obtained  in  a  retrieval  scheme  in 
degrees  of  temperature  or  percent  of  density  using  this  method  to  calculate 
transmittance  will  undoubtely  depend  on  the  number  of  sets  of  expansion  coeffi¬ 
cients. 

With  some  additional  work  two  numerical  codes  could  be  developed  to 
retrieve  atmospheric  temperature  and  pressure,  and  minor  species  concentration 
profiles  from  the  stratosphere.  These  codes  would  include  the  McMillin  method  of 
calculating  transmittance,  which  is  such  an  improvement  in  computational  speed 


that  perhaps  the  retrieval  code  could  operate  in  real  observing  time.  (Develop¬ 
ment  and  verification  of  the  codes  would  need  to  be  done  on  a  large  computer  but 
the  final  versions  could  easily  run  on  a  VAX  sized  machine.)  If  these  retrieval 
codes  were  tailored  to  the  band  pass  and  field  of  view  of  the  CIRHIS  instrument, 
routine  profiles  of  minor  and  trace  gases  in  the  stratosphere  and  on  a  gobal 


scale  could  be  obtained  from  future  CIRRIS  missions. 


APPENDIX  A 


The  geometry  of  the  limb  radiance  calculations  is  shown  in  Figure  1 .  The 
atmosphere  is  assumed  to  be  layered  and  spherically  symmetric.  Temperature, 
pressure,  and  constituent  densities  are  specified  each  kilometer  from  100  km 
down  to  the  surface.  The  horizontal  distance  from  the  tangent  point  to  some 
point  along  the  ray  path,  s,  is  related  to  altitude,  z,  as  indicated  in  Figure 
1 . 

The  "optical"  depth  is  found  by  the  expression: 

t(v,  s)  -  J  aQ  n(s)  fy(x,  y)  ds  (l) 

where  =  absorption  cross  section  at  line  center 

v  *  wave  number  at  line  center 
o 

v  «  wavenumber  (frequency) 
x  =  (v  -  vo)c0 
aD  =  doppler  half  width 

y  "  “L/aD 

*  Lorentz  half  width 
n  *  absorbing  species  density 
and  f  is  the  Voigt  line  shape  function. 

The  transmission  from  an  element  at  s  to  the  detector,  which  is  assumed  to  be  at 
s  *  <*>,  is 

Trv(s)  -  e"^w’  s)  (2) 

The  mean  transmission  over  the  wave  number  interval  v2  -  for  a  rectangular 
slit  function  is 

<Tr(s) >  -  /  2  Tr  (s)  dv/(v2  -  v, ) 

J  v, 

The  spectral  radiance  at  the  detector  is 

I(v)  -  f  *(v  -  O  B(v,  T(z))3Trv^ds 

L 


Chahine  suggested  a  technique  for  determining  temperature  or  constituent 


> 


densities  with  a  vertical  looking  instrument  which  looks  down  into  the  atmos¬ 


phere  and  scans  the  wavenumber.  The  limb  scan  problem  is  different  in  that  we 
have  a  set  of  equations  (for  a  single  line) 


W' 


gTr 

B(T(z) )  - -  ds 

as 


where  the  index  i  has  been  added  to  equation  (3)  to  refer  to  the  different 
tangent  height  measurements.  Some  of  the  constants  are  absorbed  in  Tr^.  This  is 
not  a  single  integral  equation  but  a  set  of  equations,  one  equation  for  each 
tangent  height.  Suppose  we  know  n(z)  and  wish  to  find  the  temperature  from  (4). 
The  method  of  solution  is  to  iterate  the  temperature  profile  in  the  equation 


I 


3Tr. 

B(Tm  (z))  - -  ds 

3s 


until  the  set  of  residuals 


I  -  I.” 

,m  i  i 


where  is  the  measured  intensity  at  tangent  altitudes  i  approach  a  stable  set 
of  minimum  values.  The  technique  of  iteration  is  to  satisfy  equation  (4)  at  each 
tangent  altitude  separately.  Thus  if  we  multiply  the  right  side  of  (5)  by  I^/I9 


o 


3Tr 

B(T  m  (z))  - -  ds 


and  the  temperature  is  iterated  by  the  recursive  formula 

Bd^’U))  -  B(T*  (.))  (6) 

V 

Equation  (6)  is  a  set  of  recursive  equations  giving  new  temperature  profiles 
T9*1  in  terms  of  old  profiles  T9.  We  make  the  solution  of  these  equations  unique 
by  using  the  equation  for  the  ith  tangent  altitude  to  find  the  new  temperature 


t 


e 


at  that  altitude  in  terms  Of  the  old  temperature  at  that  altitude 


B(Tm+(z  ))  -  B(t“  (z  ))  — ~ 
l  l  .  m 


where  is  the  ith  tangent  altitude.  Now  we  can  find  a  new  set  of  radiances 
Im+1  using  the  new  temperatures  (equation  5)*  The  process  continues  until  the 
residuals  become  constant  and  hopefully  small. 

To  solve  for  densities  it  is  convenient  to  integrate  (4)  by  parts  giving 

*i  *  "  Tr^(n(-*))B(T(-“)  +  j  Tr^nCz))  ds  (7) 

J  —co 

where  we  have  shown  Tr^  as  an  explicit  function  of  constituent  density  n(z). 
Assume  the  temperature  profile  is  known  or  has  been  determined  from  another 
measurement.  The  problem  is  formulated  in  much  the  same  way  as  the  temperature 
inversion  problem.  The  density  is  iterated  in  the  equations 


(j)  «  Tr(nj(-))B(T(®))  -  Tr(nj(-.  ))B(T(-»)) 


f 


Tr(n(j)(z))^  ds  (6) 
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until  the  set  of  residuals 


m  -  T  in 
ri  Ai  i 


approach  a  set  of  minima.  To  find  the  new  density  profile  in  terms  of  an  old 
profile  the  recursive  relations  are  written  as 

nn+1  ( z^)  -  aim  tx{z^  (9) 

and  the  are  found  by  some  nonlinear  technique.  One  technique  is  to  use 

Newton's  method  and  consider  the  iteration  one  tangent  altitude  at  a  time,  i.e. 

nm+1  (zi)  -  ain  n(z) 

We  start  with  -  1  and  calculate  the  corresponding  1^ ,  and  r^ .  Then  we 

2*1  2 
assume  a  second  set  of  a^,  say  a^  *  I^/I^  an<*  ca^cula^e  the  corresponding  1^ 

2  3  3 

and  r^  .  Now  in  the  linear  approximation  r^  will  be  zero  if  a^  is  given  by 


m 


m-2  m-1 

m-2  m-2  ®i _ ~al 

ai  "  °i  ~  ri  m-2  m-1 

ri  "ri 

This  formula  can  be  used  for  further  iterations  using  for  r^m  ^  the  minimum  of 
all  previous  r^°  and  the  corresponding  a^m  for  The  process  terminates 

when  the  set  of  r^m  are  minimum.  Having  found  a  set  of  a^  we  now  return  to 
equation  (9)  and  find  a  new  n(z)  and  repeat  until  convergence  is  obtained  in 
the  sense  the  3^  all  approach  one. 


